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We present here an overview of the history, applications and important properties of a function 
which we refer to as the Levy integral. For certain values of its characteristic parameter the Levy 
integral defines the symmetric Levy stable probability density function. As we discuss however the 
Levy integral has applications to a number of other fields besides probability, including random 
matrix theory, number theory and asymptotics beyond all orders. We exhibit a direct relation¬ 
ship between the Levy integral and a number theoretic series which we refer to as the generalised 
Euler-Jacobi series. The complete asymptotic expansions for all natural values of its parameter are 
presented, and in particular it is pointed out that the intricate exponentially small series become 
dominant for certain parameter values. 

PACS numbers: 


Dedicated to the memory of Elliott Montroll 


I. INTRODUCTION 

The objective of this article is to pay homage to an integral which has occupied the minds of a host of great 
mathematicians including Cauchy, Levy, Hardy, Littlewood, Ramanujan, and Polya and which continues to arise in 
modern contexts such asymptotics beyond all orders and random matrix theory. This function, Fa{z), which we refer 
to as the Levy integral is defined via 

poo 

Fa{z) = / exp{—t°') cos{zt) dt a > 0. (1) 

Jo 

For a < 1 the integral is defined only for real z. 

The number of seemingly disparate topics in which Fa{z) appears in a fundamental way is remarkable. The reader 
familiar with the theory of probability will immediately recognise the function (l/7r)Fh(z) as defining the symmetric 
Levy stable probability density function when a <2. For a = 1,2 the function {l/TT)Fa{z) is identified respectively 
with the Cauchy and Caussian probability densities. However, Fo,(z) also has an intimate connection to Waring’s 
problem in number theory, is of great interest in asymptotics beyond all orders, and appears in certain problems 
arising in random matrix theory and wave phenomena. 

The analytic properties of Fa{z) are of long standing interest. Bernstein in 1919 investigated the zeros of Fh(z), 
showing that it possesses an infinite number of real zeros when a = 4,6,8,.... He appears to have been the first to 
realise that Fa{z) might become negative for some values of z when a > 2. In such cases {l/Tr)Fa{z) does not dehne 
a probability density function. In 1923 Polya ||, ^ extended the work of Bernstein by systematically investigating 
the zeros of iC,(z) and showed that when a = 4,6,8,..., Fa{z) has an infinite number of real zeros but no complex 
zeros, whereas if a > 1 is not even integer then Fa{z) has an infinite number of complex zeros and a finite number, 
not less than 2[a/2], of real zeros. Levy [Q in 1923 proved the non-negativity of Fa{z) for all real z when 0 < a < 2, 
and was responsible for the development of the theory of stable distributions [^, and so we refer to the function 
Fc,{z) as the Levy int^ral for all positive a regardless of whether it defines a probability density function or not. It 
has been observed that for all positive rational a the function Fq(z) can be expressed as a sum of generalised 

hypergeometric functions, and in fact Fa (z) is itself a specific type of Fox H-function . There are however only a 
handful of specific values of a for which an identihcation with more familiar special functions is known 

The present authors 0 have recently constructed the complete asymptotic expansions of Fa{z) for all natural a, 
and in particular have demonstrated that there exist very intricate series of exponentially small terms lying beyond all 
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orders of the asymptotic power series. Interestingly, when a is an even integer the asymptotic power series vanishes 
and these exponentially small series are liberated from their subdominance. This provides a powerful case study 
for the newly emerging field of asymptotics beyond all orders |l^, the discussion of which we shall return to after 
discussing how Fa{z) arises in the fields of probability, random matrix theory, diffraction integrals, and number theory. 


II. LEVY STABLE PROBABILITY DENSITY FUNCTIONS 

The theory of stable distributions is a branch of probability theory of fundamental importance. At a time when 
it was fashionable to find the most general conditions possible in order for the central limit theorem to be satisfied, 
Paul Levy was interested in finding elegant examples of when it is not. Indeed Levy ’s stable probability density 
functions have divergent moments for a ^ 2 , and the definition of stability is manifestly one of self similarity rather 
than universality. Levy defines a probability density function p{x) to be stable if it satisfies the following property 
Si a- Let Xi and X 2 be independent random variables with the same probability density function p{x), then 

ciAi + C 2 A 2 = cA ( 2 ) 


for all positive ci and C 2 , where X has the same probability density function as Xi and A 2 . The constant c is 
determined by the constants ci and C 2 . The symmetric solutions of this condition have characteristic functions of the 
form exp(—|g|“) Q, which is the Fourier transform of {l/'!r)Fa{z). 

The function Fa{z) was studied by Cauchy in the early 1850s in relation to the theory of errors, and it is 
evident that he was aware that these functions satisfy Levy stability. Cauchy appears unaware however that Fa{z) 
may become negative for some values of z when a > 2 and so fail the non negativity requirement of a probability 
density function. 

The divergence of the moments of F^^z) for a not an even natural can immediately be seen from the non-analyticity 
of exp(—|g| ), since if the characteristic function of a probability density function is analytic then it possesses a power 
series whose coefficients are proportional to the moments of that probability density function. We can summarise[^ 
the behaviour of the moments of (l/ 7 r)Fh(z) as a function of a as 


< z 


2m 


> = 



a = 2fG {2,4,6,...; 

J'- 

+00, 

a < 2m and a G U“ 

— OO, 

a < 2m and a G U“ 

0 , 

otherwise. 


(4/c,4fc + 2) 

(4fc + 2,4fc + 4) 


( 3 ) 


The stable probability density functions have been experiencing a tsunami of popularity in a wide variety of applica¬ 
tions in the past decade, from physics to biology to economics |l^. The reason for their utility stemming precisely 
from the divergence of the moments for a 2. This implies that for a < 2 no characteristic scale exists and random 
walks constructed using (l/ 7 r)FQ,(z) as the step length probability density function display statistically self similar 
trajectories, and as stressed by Mandelbrot the fractal dimension of these trajectories is |^, |^. Thus Fa{z) 
plays a fundamental role in the increasingly popular pursuit of modeling phenomena which display structure on many 
scales. 


III. RANDOM MATRIX THEORY AND THE CUSP DIFFRACTION CATASTROPHE 


It has been noted Q that ^ 4 ( 2 /) is proportional to the special case of the Pearcey integral P{x,y) in which its first 
argument is set to zero. Specifically, 


P(0,y) = 2e*-/8F4(y). 


( 4 ) 


The Pearcey integral is the canonical form of the oscillatory integral describing the cusp diffraction catastrophe 
^ and is defined by 


/ OO 

exp[i(u'^ -I- Xu^ + Y u)] du, 

-00 

which can be expressed in an alternative more general form via analytic continuation 

poo 

P{x,y) = / eyip{—t‘^ — xt^)cos{yt)dt. 

Jo 


( 5 ) 


( 6 ) 
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This integral was investigated by Pearcey in his investigation of the electromagnetic field near a cusp [^if and is utilised 
in many short wavelength problems such as wave propagation and optical diffraction 11, 00- The function P{x, y) 
was known before Pearcey however, and was discussed by Lord Rayleigh in 1879 in connection with the diffraction of 
light Brillouin, as early as 1916, obtained an asymptotic expansion of P{x,y) using steepest descents [^. As 
discussed in ||^, P{0,y) displays dominant exponentially small oscillatory behaviour. As we shall discuss later, this 
behaviour is shared by Fa{z) for all even natural a larger than 2. 

Brezin and Hikami have recently demonstrated that Fi{z) appears in the solution of a certain problem in 
random matrix theory. They consider an N x N Hamiltonian matrix, F[ = FIq + P, which is the sum of a given 
nonrandom Hermitian matrix, iJo, and a random Gaussian Hermitian matrix V. The addition of FIq breaks the 
unitary invariance of the probability measure that defines the Gaussian unitary ensemble, 

PiH) OC OC e-W2)TriH--2HoH) ( 7 ) 


They investigate the eigenvalue density of when the eigenvalues of the nonrandom matrix, iJo, are chosen to be 
±a, each value being N/2 times degenerate. In this case when a > 1 a gap appears in the eigenvalue spectrum of F[ 
which closes as a ^ 1. 

An important quantity in random matrix theory is the n-point correlation function, A 2 ,..., A„), given by 

R„(Ai, A 2 ,..., A„) = > (8) 


where \j is an eigenvalue of FI. 

The presence of nonzero Hq implies that the usual orthogonal polynomial methods are no longer available for the 
determination of the n-point correlation function. Brezin and Hikami show however that the function i?„(Ai, A 2 ,..., A„) 
is again expressible as the determinant of an n x n matrix whose elements are given by a kernel K^{Xi^ Xj), just as 
in standard random matrix theory. The following contour integral expression for the kernel is arrived at 


/rN(Ai,A.) = (-ijww r^i 

J -00 27r J\u\>a 27rz \ 


du f 


(9) 


1 


,-{N/2)u^-{N/2)t^-NitXi+Nu\2 


U — it 


For the critical value of a = 1 Brezin and Hikami then find that in the scaling limit in which N is large and A^^/^A 
is finite, the kernel becomes 


K{x,y) = 


<j)'{x)i)'[y) - ^"{x)ip{y) - $(x)ip"{y) 


where 


and 


x-y 


^{x) = —F4(-\/2a:), 

TT 




( 10 ) 


( 11 ) 


( 12 ) 


They then go on to show that this behaviour is universal in the sense that K{x,y) is independent of the eigenvalue 
distribution of the non random Hermitian matrix Hq provided that a gap closes at the origin in the eigenvalue density 
of the total matrix H. 

Brezin and Hikami then discuss a situation in which one considers a random Hermitian matrix coupled to a non 
random complex matrix, Hq, rather than an Hermitian matrix. In this case the probability measure of the total matrix 
becomes complex and so loses its probabilistic interpretation. By appropriately selecting the eigenvalues of this non 
random matrix Brezin and Hikami discuss a multicritical generalisation of the closing gap phenomena discussed above, 
the kernel of which is now expressible in an analogous manner in terms of F 2 k{z), with each natural k larger than 2 
corresponding to a particular choice of the eigenvalues of Hq . As we mentioned in the introduction, and will discuss 
in more detail presently, it is precisely this subclass of the functions Fa{z) that display dominant transcendental 
asymptotics. 
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IV. THE LEVY INTEGRAL, THE GENERALISED EULER-JACOBI SERIES AND WARING’S 

PROBLEM. 

Waring’s conjecture asserts that every natural number can be expressed as the sum of s integers all to the kth 
power. Let the function g{k) denote the least value of s such that the proposition is true for all natural numbers, 
and G(fc) the least value of s such that the proposition is true for all but a finite number of natural numbers. 
Waring’s problem naturally decomposes into two parts, that of proving the existence of g{k) for all fc, and that of 
the determination of the explicit structure of g(fc) and G{k). Hilbert|M proved the existence of g{k) in 1909, and 
although much progress took place last century the investigation of g{k) and G{k) continues to this day. 

A fundamental arithmetic function that arises in the investigation of Waring’s problem is which for 

k even is the number of representations of n G N in the form 

n = a!l + a\+ all ai^'L, (13) 

with representations which differ only in the order of the reckoned as distinct. The generating function of the 
arithmetic function rk^sin-) which was investigated by Hardy and Littlewood in relation to Waring’s problem, 
takes the elegant form 


Y,rkAn)e-‘^’^ = [2Sk{a)-ir (14) 

n—0 

where the function S'fe(a) is special case of a function which we refer to as the generalised Euler-Jacobi series, defined 
by 

oo 

S'!(a) = ^ p,gGN 3fJ(o)>0. (15) 

n—0 

The series Sp{a) was studied in detail in |^. Application of Euler’s summation formula |^, results in the 

following inversion, which is the generalisation of the familiar transformation identity for the Jacobi theta function 
93(a) = 252(0) - 1 . 


Se (a) = 


T (q/p + 1) 

nllP 


■E 

n=l 




2mT 


p/i I „«/p 


(16) 


Hence we see that there exists a startlingly direct relationship between the number theoretic Euler-Jacobi series 
Sp/q{a), and the Levy integral Fpiq(z). For this reason, Hardy and Littlewood devoted an entire section of [p^ to 
the investigation of Fa{z), where it is noted that Fi(z) can be written as a sum of hypergeometric functions and that 
it displays dominant exponential asymptotics, two observations which we show are indeed true for all even natural 
a. Without the modern theory of the generalised hypergeometric functions at their disposal, as appears in [ p^ for 
instance, an investigation along the hypergeometric route however is dismissed by Hardy and Littlewood as containing 
too many formal complications. In 0 the hypergeometric representations are given for Fp/q{z) and hence Spfq(a), 
and the complete asymptotic expansions are investigated for various cases of Ramanujan was interested in the 
series Sp/q{a) and derived a small a asymptotic expansion for it as discussed by Berndt |30|] . This expansion consists 
only of the component power series however which vanishes identically when a is an even natural, all exponentially 
small terms having been neglected. Berndt uses the Mellin transform technique rather than utilising the inversion 
formula (16) and hypergeometric representation of Fa{z) to recover Ramanujan’s result. 


V. THE EXPONENTIAL ASYMPTOTICS OF THE LEVY INTEGRAL 

The question of the asymptotics of Fa{z) appears to have been first investigated thoroughly by Burwell |3|], who 
utilised the technique of steepest descents to obtain the first terms in the dominant asymptotic expansion but was not 
concerned with obtaining the exponentially small subdominant terms. By expressing Fo,(z) as a sum of generalised 
hypergeometric functions the complete large 2 : asymptotic expansions of Fa(z) were calculated in [|| for a = 3,4, 5, 6, 7, 
which then yielded the small a asymptotic expansions of Sp/q(a). A more general problem has recently been solved 
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by the present authors. The d dimensional isotropic Levy stable probability density function possesses the following 
integral representation, where J(d/ 2 )-iiQ'f) is a Bessel function, 

p» = tAf / d'ge—'iexp(-lqr) (17) 

oo 

= (2^^)(d/2) (18) 

0 

When d — 1 this clearly reduces to {l/7T)Fa{z). It is shown in 0 that (r) can be written in terms of generalised 
hypergeometric functions for all positive rational a. These hypergeometric representations are then used to construct 
the complete asymptotic expansions for ([l8| ) for all natural a. Utilising the useful notation, /(±z)=*= = f{z) ± f{—z) 

for some function /, the Levy integral for rational a can be identified with a sum of generalised hypergeometric 

functions as follows 




X 



(2^) 


(p-g-2) 

2 


p-1 


1=0 


n(P-Lr[(M 



21 



2i + l I 1 21 + 1 I 2 2^4-1 I q-1 

’ p ■'"g’ p p 'T g 

2 Z + 2 2 i + 3 21+p 

p • p V, p 






+ 


(19) 


where p, g are natural numbers. The complete asymptotic expansions of Fa{z) and S'a(a) for all natural a are a 
straight forward corollary of the asymptotics for p(( (r) given in 0, and we present the expansion for Fa{z) below. In 
this way, this article can be considered as an update of the work presented in 0 , in that we now display the structure 
of the asymptotics for all a G Z> 3 , rather than discussing specific cases. We remark that the complete asymptotic 
expansion of Pai^) presented in [ |l0| is only extremely weakly dependent on the dimension d, none of the overall 
structure changing as one moves from d = 1 to d > 1. Thus all of the interesting asymptotic behaviour possessed by 
p^(r) is retained by Fa{z) and in this sense the asymptotic behaviour ofp^(r) for natural a is universal with regard 
to the spatial dimension d. 


Fa ( 2 ) - 


X 


X 


X 


— (—)"*+^ r(ma) , /rmra\ 


-U 


m —1 


r(m) 
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sin 




2 (a - 1) (az«-2)l/2(a-l) Z! 2 ( z ) 

^ ^ ^ ' m —0 

[a/2] - [a/4]-1 

exp I —(a — 1) sin 


Q;\"*“/(a-l) 
z - 


k=0 

COS I (a — 1) cos 


(4fc + 1)74 


[2(0-1) J 


Z \ a/(a-l) 


o 


{Ak + 1)74 

2(a-l), 

^(a—4[a/4]),2 


2(a- 1) (az“-2)i/2(a-i) 


(4/4 + 0)77474 ^ (4/4 + 2 — 0)74 


771=0 


2(0-1) 4(0-1) 

'o\"*“/(“-i) 

. z - 


i-y 


-(„-i)(i) 




exp 

TT / X TT 

- < arg(z) < 
a a 


00, 


CX G ^>3- 


( 20 ) 


Here [x] refers to the integer part of the real number x, and 5i^a-i[a/A]),2 is the Kronecker delta. The Nk are defined 
by the following recursion relation 
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2 (a-l) 2(a-l)-s 

s—1 r=0 


2(a- 1) 


—-— r + s - fc) - a - I + —-- , 

2{a-l) 


(r + s — fc) — (a — I) — 


()«+.[2(a-l)]2 


rl[2{a — I) — s — r\\a^°‘~^k 


(«- 2 ) 
4(a - I) 




iV. = 0 


s < 0, 


iVn = I. 


( 21 ) 


where {x)a denotes the Pochhammer symbol. 

We note that for the case of a = 4, is in agreement with the corresponding result for the Pearcey integral given 
in We remark also that in the process of extending Polya’s work on the zeros of Fa{z), Senouf |3^ has obtained 
the fc = 0 portion of the complete asymptotic expansion for even a. 

The small a asymptotics of the series Sa{a) is now given simply by utilising its relationship with the Levy integral 
&■ Since the cases a = I, 2 are simply the Cauchy and Gaussian probability density functions, up to a normalisation 
factor of TT, whose asymptotic behaviour is trivial to obtain, ( p0[ ) gives us a complete understanding of the asymptotic 
behaviour of Fa{z) for all natural a. It is clear from the form of the power series component of the full asymptotic 
expansion, given by the first series of (20), that when a is even every term vanishes identically and so the exponentially 


small terms become dominant. We note that when a is an even integer the hypergeometric representation of Fa{z) 
simplihes in such a way that only hypergeometric functions with no numerator parameters appear |p^ , and it is known, 
see e.g. [^, that such hypergeometric functions always display purely transcendental asymptotics. Thus the function 
Fa{z) displays qualitatively different behaviour depending on whether the parameter a is even or odd. Also, the 
presence of [a/4] in the summation limits and the Kronecker delta show that the structure of the asymptotics of 
Fa{z) depends sensitively not only on the parity of a, but on the arithmetic residue of a modulo 4. The mechanism 
by which this dependence arises in the construction of (|2^) is discussed in [^ . An interesting property of (|^) is that 
the number of exponentially decaying oscillatory series increases as a increases, so that the asymptotics becomes ever 
more intricate as a becomes large. This phenomena was observed in |^, for the specific cases a = 3,4, 5,6, 7, and 
predicted to be true for general a. Our result (p0|) now confirms this. 


VI. CONCLUDING REMARKS 

The theory of stable distributions is a fundamental branch of probability, and consequently the importance of the 
Levy integral for a less than two is widely recognised. What is equally as intriguing as its probabilistic applications 
is the fact that the Levy integral also appears to be of importance when a is larger than 2 and it no longer defines a 
probability density function. Its delicate asymptotics when a is an even natural number is qualitatively different to 
anything exhibited when a < 2. From examples such as its appearance in random matrix theory, and its relationship 
with the Pearcey integral, one might well expect that the future growth of interest in the Levy integral might lie 
not in probability, in which it is well established, but in problems in which its exponential asymptotics are utilised. 
In particular it appears that the subclass of Fa{z) where a is an even natural number, which exhibits dominant 
exponential asymptotics, will be of interest. 
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